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We study the coherent cooperative phenomena of the system composed of two interacting atomic 
ensembles in the thermodynamic limit. Remarkably, the system exhibits the Dicke-like quantum 
phase transition and entanglement behavior although the governing Hamiltonian is fundamentally 
different from the spin-boson Dicke Hamiltonian, offering the opportunity for investigating collective 
matter-light dynamics with pure matter waves. The model can be realized with two Bose-Einstein 
condensates or atomic ensembles trapped in two optical cavities coupled to each other. The inter- 
action between the two separate samples is induced by virtual photon exchange. 
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I. INTRODUCTION 

Entanglement, arising from nonlocal superposition 
state of two or more quantum systems, is one of the most 
striking features of quantum mechanics and plays a key 
role in the test of local hidden variable theories [1,2]. 
Over the past two decades entanglement is regarded as a 
key resource to implement quantum information tasks, 
such as quantum cryptography [3], computer [4], and 
teleportation [5]. On the other hand, entanglement is 
closely related to various collective quantum phenomena 
and plays a central role in studying quantum many-body 
systems. Typical examples are entangled ground states 
used to clarify superconductivity [6] and fractional quan- 
tum Hall effect [7]. Ground state entanglement, which is 
responsible for long-range correlations at zero tempera- 
ture [8], may shed light upon quantum phase transition 
[8]. The connection between entanglement and quan- 
tum phase transition has been extensively explored in 
the quantum system made up of N spins 1/2 on a one- 
dimensional lattice [9-11] or on a simplex [12,13]. 

Another important example of correlated many- 
particle systems is the Dicke model [14], which describes 
the interaction of an ensemble of N two-level systems 
with a single-mode quantized field. When the atom-field 
coupling is strong enough such a model exhibits a super- 
radiant quantum phase transition [15-17]. Lambert et al. 
have investigated the critical behavior of the atom-field 
entanglement and pairwise entanglement between atoms 
in the thermodynamics limit [18]. The Dicke model has 
also been a fertile ground for implementation of quan- 
tum information as atomic ensembles with long-lived 
electronic states are ideal for storing and processing lo- 
cal quantum information via interaction with light fields. 
High-fidelity quantum operations and entanglement may 
be achieved beyond the strong coupling between a single 
atom and a single photon due to the collective enhance- 
ment of interaction with the field mode [19]. 

The Dicke quantum phase transition occurs under the 
condition that the coupling strength is on the order of 



the energy separation between two involved levels and 
thus the counter-rotating terms significantly affect the 
dynamics, which is difficult to achieve in typical cavity 
QED systems since the atomic transition frequencies usu- 
ally exeed the atom-cavity coupling strength by many 
orders of magnitude. Dimer et al. have propsed an effec- 
tive Dicke model based on balanced Raman transitions 
between two ground atomic states of an atomic ensem- 
ble induced by a cavity mode and a pair of laser fields 
[20], which significantly lowers the energy difference be- 
tween the two involved levels. It has been shown that the 
Dicke model can also been realized in a laser-driven Bosc- 
Einstcin condensate coupled to an optical cavity and the 
onset of self-organization corresponds to the Dicke quan- 
tum phase transition [21-23]. However, these systems are 
subject to cavity loss, which would significantly alter the 
coherent dynamics of the Hamiltonian and deteriorate 
the matter-light entanglement. Though the superradiant 
quantum phase transition has been observed in Ref . [22] , 
the ground state entanglement has not been verified. 

In this paper we investigate the collective dynamics of 
the system that consists of two coupled atomic ensem- 
bles and find two main results. Firstly, we show that 
the coupled spin system can display the quantum phase 
transition and vacuum entanglement properties of the 
physically distinct spin-boson Dicke model in the ther- 
modynamic limit, providing an access to the intriguing 
critical entanglement of the Dicke model without involv- 
ing the dynamics of the quantized field. Secondly, as an 
example for the physical implementation of this model we 
demonstrate that two Bose-Einstein condensates (BECs) 
or atomic ensembles trapped in two coupled optical cav- 
ities can form the effective spin Hamiltonian. With the 
assistance of suitable chosen external fields, the virtual 
excitation of the cavity modes mediates effective coupling 
between the two spin ensembles. Besides fundamental 
interest, the system offers the possibility to produce two- 
mode squeezed vacuum states between two spatially sep- 
arated macroscopic systems, which are of crucial impor- 
tance for quantum communication [24] and nonlocality 
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test [25] with continuous variables. Such states are also 
useful for exploring the boundary between quantum and 
classical worlds and understanding the decoherence effect 
in quantum information. 

The paper is organized as follows. In Sec. 2, we study 
the Hamiltonian dynamics of two coupled atomic ensem- 
bles in the thermodynamic limit, and show that the sys- 
tem can exhibit the quantum phase transition and criti- 
cal entanglement behavior of the Dicke model. In Sec. 3, 
we propose an experimental realization of the system by 
coupling the motional degrees of freedom of two BECs 
trapped in two coupled optical cavities. The model can 
also be realized in the electronic degrees of freedom of 
two atomic ensembles via virtual photon exchange. A 
summary appears in Sec. 4. 



II. QUANTUM PHASE TRANSITION AND 
VACUUM ENTANGLEMENT 



Let us start by considering the system involving two 
coupled atomic ensembles, with the ith (i = 1, 2) ensem- 
ble consisting of two-level atoms of energy splitting 
uji . The Hamiltonian of the system is given by 
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Under the critical point A c = the solution is 

a = (3 = 0, which corresponds to the normal phase. 
In this case the effective Hamiltonian = u>ib\b\ + 

u 2 b\b 2 + X(b\ + W){b\ + b 2 ) is mathematically equiva- 
lent to the spin-boson Dicke Hamiltonian in the normal 
phase [15]. The ground state in this phase is a two-mode 
squeezed vacuum state, which has a definite partiy and 
is incoherent, i.e., ((Jj + + J^)) = 0. 

Above the critical point, there exist two physically 
sensible solutions {a — = —\/N 2 j3o} an d 

{a = —^/Nxa 0l (3 — a/A^A)}, where a — {^[1 — 

{i[l- 
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and r = N 2 /N\. The corresponding effective Hamilto- 
nian is 
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H = uji J\, z +oj 2 J 2 ,. 



-^= { Jt + J^Ji + J 2 ), (1) 

where {Ji, z ,J^} are collective atomic operators for the 
ith ensemble, satisfying the angular momentum commu- 
tation relations [J* , J^} — 2Ji.z and [J^, Ji, z ] — -fjf- 
The Hamiltonian commutes with the partity operator 
e *7r(./i,z+./2,*) Using the Holstein-Primakoff representa- 
tion, we express the operators {Ji. Zl J^ 1 } in terms of the 
annihilation and creation operators bi and b\ of a bosonic 

mode via J+ = b\\J - b\bi, J[~ = \J~Nj~- b\bibi, and 

Ji,z = b\bi — Ni/2. Then we obtain the two-mode bosonic 
Hamiltonian 



H = u>ib\bi + uj 2 b\b 2 



>I yj^-blh + \j N x - b\bib x 
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Perform 



the 



transformation 
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D+(a)D+ ((3)HD 1 (a)D 2 ((3), where Di(oi) = e' 
and D 2 (P) = e l3b i^ l3 * hi are displacement operators. In 
the thermodynamic limit JVj -) oo we can approximate 
the square root terms up to l/Ni. The resulting Hamil- 
tonian is expanded up to the second-order in the boson 
operators. 

In order for the coefficients of the linear terms in 
the displaced Hamiltonian to be zero, the displacement 
amounts a and /3 should satisfy 
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The Hamiltonian H ^ is diagonal in terms of two normal 
bosonic modes with the frequencies 



w, = 



K X K 2 + K 3 K 4 



± 



(K X K 2 - K 3 K 4 ) 2 



-,1/2 



+ 4K 2 K 4 Kl 



■ (6) 



In the regime near the critical point, the energy gap to the 

1 /2 

first excited state vanishes as cj_ oc |A-A C | 1/Z . Above 
the threshold, in addition to two-mode squeezing both 
modes exhibit single-mode squeezing since the coefficient 
of b\ 2 + b\ docs not vanish for each mode. However, 
the single-mode squeezing does not affect the entangle- 
ment between the two collective atomic modes. Below 
the threshold, the single-mode squeezing vanishes and 
thus the threshold is not affected. 
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Above the threshould, there are two degenerate ground 
states, which corresponds to the breaking of the parity 
symmetry with the atomic polarizations ((Ji + Ji)) = 
2ay/Nx - a 2 and (( j£ + J 2 ~)) = 2(3y / N 2 - (3 2 acquiring 
macroscopic populations. In this phase both the coher- 
ent excitation and squeezing contribute to the excited 
atomic numbers. This is analogous to the superradiant 
phase of the spin-boson Dicke model with the field mode 
replaced by another atomic mode. The excitation num- 
ber of each ensemble as a function of A/A c is ploted in 
Fig. 1. The solid line represents the incoherent excita- 
tion number due to squeezing, while the dashed line rep- 
resents the scaled coherent excitation number. Near the 
critical point the incoherent excitation number of each 

atomic ensemble due to squeezing is n inc oo |A — AcP 1 ^ 4 , 
which diverges with the same exponent 1/4 as the cor- 

— 1/2 

relation length £ = lu_ . The entanglement between 
the two atomic ensembles can be determined by the von 
Neumann entropy S = — trplog 2 p, where p is the re- 
duced density operator of one atomic ensemble. In the 
superexcitation phase, this entropy is given by 



S = (k + 1/2) log 2 (fc + 1/2) - (k - 1/2) log 2 (fc - 



where k = 
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1/2 

and 



6=\ 'AYCt&n^KsyfK^Kl/iK^ - K 3 K 4 )}. The entan- 
glement as a function of A/A c can be seen in Fig. 2. At 
the critical point the entanglement diverges logarithmi- 
cally also with the exponent 1/4 as 5oo log 2 |A — A c | -1 ^ 4 , 
which is analogous to the critical behavior of the atom- 
field entanglement in the Dicke model [18]. In the above 
we have assumed that A > 0. It should be noted that the 
system exhibits the same critical behavior and entangle- 
ment for A < 0. 



III. PHYSICAL REALIZATION 

We note that the Hamiltonian (1) can be realized in 
the quantum motions of two BECs or atomic ensem- 
bles trapped in two coupled single-mode cavities along 
the x-axis. The resonant coupling between the two 
cavity modes is given by the interaction Hamiltonian 
H c = v{a\a2 + 0,10,2), where a\ and dj are the creation 
and annihilation operators for the ith cavity mode, and 
v is the intercavity hopping strength. Such a coupling 
can be mediated by overlap of their evanescent fields or 
by an optical fiber [26]. 

We first consider the cavity-BEC system. Suppose that 
the ith BEC is composed of Ni two-level atoms, each of 
which is coupled to the ith cavity mode with the max- 
imum coupling strength gi and driven by a pump laser 
field along the j/-axis with the maximum Rabi frequency 
fli. The pump frequency oo p , close to the cavity mode 
frequency ui c , is highly detuned from the atomic transi- 
tion frequency uj a . Then the atomic excited level can be 



adiabatically eliminated and the atoms coherently scat- 
ter light between the pump field and the cavity mode, 
which induce two balanced Raman channels between the 
atomic zero-momentum state and the symmetric super- 
position of states with the momentum of a photon, Hk, 
along the x and y directions. Hence each atomic field 
can be described in a Hilbert space spanned by two 
Fourier-modes cq.i and ci^, with cj ^0,1 + c\ ^Ci,, = Ni 
being a constant of motion [21,22]. The coupling be- 
tween the two modes can be described by the angular 
momentum operators: J 4 + = c\ ^0,1, J~ = c' oi Ci t i, and 

Ji,z — 2"( c ii c M ~ 4 t c o,i)- I R the frame rotating with 
the pump field frequency co p the Hamiltonian for the ith 
BEC-cavity system is 

Hi = 8 c ^a\a i +LUQj i ^~^{ai+aX)(J++J^)-Uia\a i c\ l ci il 

'(8) 

S c ,i = Uc-u p -Nigj /(2A a ), uj = hk 2 /m, rji = gjOj/A a , 



3<??/(4A a ), and A a 



U}„ 



bj c . Under the con- 



dition N\ — N 2 and g\ = 52 we can choose the pump 
frequency appropriately so that S c ,i = 0. Introducing the 
new bosonic modes di.2 = ^-(ai ± a 2 ), we can rewrite 
the Hamiltonian of the total system as 

H t = vd\d\ - vd\d 2 + ujq J\,z + ^n-h.z 



(d 1 +d\ +d 2 +d\)(J+ + J 1 
(di + d\ - d 2 - d\){J$ + J2 



m 

2V2 
m 

2V2 

- y (<4 rf l + d 2 d 2 + d l d 2 + 4 d l) c l,l c l,l 

--^{d\d\ + d\d 2 - d\d 2 - rf 2 ^i) c i,2 c i,2-(9) 



Under the condition v ^> \fN\r\i, ujq, m, 5 c ,i , where 

rij being the mean motional excitation number in the ith 
BEC, one can adiabatically eliminate the field modes to 
obtain the coupling between the two condensates. So the 
effective Hamiltonian is given by 



H e = 



vd\d\ — vd\d2 + ojoJi.z + ^0 J2, 

+ ^_(J+ + Ji)(J+ + J2) 



+ 



VKN2~ 
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p(4^i - 4 rf 2)(«icl il ci > i - u 2 cl i2 ci i2 ) 2 



4;/ 
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where A = - y/N\ N 2 mm / fa) ■ The Hamilt onian de- 
scribes a four-photon process which is induced by virtual 
excitation of the atomic electronic states and field modes. 
We note that there is no coupling between atoms belong- 
ing to the same ensembles since the detunings of the two 
nonlocal field modes d\ and d 2 are opposite, which leads 
to opposite contributions to the coupling. On the other 
hand, these two nonlocal modes equally contribute to the 
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couplings between atoms belonging to different ensem- 
bles because the product of the two Raman transition 
coefficients associated with mode d\ is also opposite to 
that associated with cf 2 . When the two cavity modes 
are both initially in the vacuum state, the two bosonic 
modes d\ and di will approximately remain in the vac- 
uum state during the evolution since their frequencies 
are highly detuned from the pump frequency due to the 
strong coupling between the two cavities. In this case 
the effective Hamiltonian reduces to Eq. (1), with the 
effective coupling strength A being controllable by the 
Rabi frequencies or detunings of the pump fields. The 
quantum phase transition corresponds to the simultane- 
ous self-organization of the two condensates. In a realis- 
tic experiment, the system is a driven and damping one, 
which will realize a steady state governed by energy flow 
from the pump fields into the cavity fields, rather than 
a true ground state of the Hamiltonian, similar to that 
studied in Rcf. [22]. 

The Hamiltonian (1) can also be realized in the 
electronic degrees of freedom of two atomic ensembles 
trapped in two coupled cavities. The cavity mode, to- 
gether with two external fields, can induce balanced off- 
resonant Raman transitions between two ground states of 
each atomic ensemble [20,27]. With appropriate choice of 
the parameters of the external fields, the field modes can 
be adiabatically eleminatcd and the two atomic ensem- 
bles are coupled via virtual photon exchange. Due to the 
stability of the atomic ground states the vacuum entan- 
glement between the two atomic ensembles should have 



a long coherence lifetime and can be readily transferred 
to light fields. 



IV. SUMMARY 

In conclusion, we have investigated theoretically the 
ground state properties of the model involving two cou- 
pled spin ensembles in the thermal limit, showing that 
the model displays the quantum phase transition and vac- 
uum entanglement described by the Dicke model despite 
the fundamental disctinction between these two models. 
The model can be realized in the motions of two BECs 
or in the internal states of two atomic ensembles in two 
coupled cavities. The coupling strength between the two 
spin ensembles can be tuned via the parameters of exter- 
nal fields, making the system a promising simulator for 
this model. The entanglement within each atomic en- 
semble and its connection with the quantum phase tran- 
sition will be further investigated. Another interesting 
problem is how the interaction between atoms belonging 
to the same ensembles affects the critical behavior and 
entanglement. 
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FIG. 1: (Color online) The mean excitation number as a func- 
tion of A/A c . The solid line represents the incoherent excita- 
tion number due to the squeezing, while the dashed line rep- 
resents the scaled coherent excitation number, given by the 
coherent excitation number devided by N. The parameters 
are wi = tJ2 = w, A c = uj/2, and Ni — N2 — N. 



FIG. 2: (Color online) The entanglement entropy between the 
two atomic ensembles as a function of A/A c . The parameters 
are wi = 0J2 = u), A c = w/2, and Ni — N2 = N. 



